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Abstract. We explore the systematics of the density dependence of nuclear 
matter symmetry energy in the ambit of microscopic calculations with various 
energy density functionals, and find that the symmetry energy from subsat¬ 
uration density to supra-saturation density can be well determined by three 
characteristic parameters of the symmetry energy at saturation density po, i-e., 
the magnitude Esym(po), the density slope L and the density curvature Ksym- 
This finding opens a new window to constrain the supra-saturation density be¬ 
havior of the symmetry energy from its (sub-)saturation density behavior. In 
particular, we obtain L = 46.7 ± 12.8 MeV and Ksym = —166.9 ± 168.3 MeV as 
well as Esym(2po) ~ 40.2± 12.8 MeV and L(2po) ^ 8.9± 108.7 MeV based on the 
present knowledge of Esym(po) = 32.5 ± 0.5 MeV, Esym(pc) = 26.65 ± 0.2 MeV 
and L{pc) = 46.0 ±4.5 MeV at pc = 0.11 fm“^ extracted from nuclear mass and 
the neutron skin thickness of Sn isotopes. Our results indicate that the sym¬ 
metry energy cannot be stiffer than a linear density dependence. In addition, 
we also discuss the quark matter symmetry energy since the deconfined quarks 
could be the right degree of freedom in dense matter at high baryon densities. 


1 Introduction 


The nuclear matter symmetry energy, which essentially characterizes the isospin dependent 
part of the equation of state (EOS) of asymmetric nuclear matter, is important for under¬ 
standing many questions in nuclear physics and astrophysics, including the nuclear effective 
interactions in asymmetric nuclear matter, the structure and stability of exotic nuclei, the re¬ 
action dynamics induced by rare isoto^s, the nature and evolution of neutron stars, and the 
mechanism of supernova explosion [l|-l2l • The symmetry energy also plays an important role 
in some interesting issues of new physics beyond the standard model (ifll. During the last 
decade, a lot of experimental, observational and theoretical efforts have been devoted to con¬ 
straining the density dependence of the symmetry energy 13 in. While significant progress 
has been made in determining the density behavior of the symmetry energy around satura¬ 
tion density po (^ 0.16 fm“^), its supra-saturation density behavior is still poorly known and 
remains the most uncertain property of isospin asymmetric nuclear matter. Theoretically, 
many experimental and observational probes have been proposed to extract information on 
the supra-saturation density behavior of the symmetry energy [H. In terrestrial laboratories. 
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heavy-ion collisions provide the only way to explore the snpra-satnration density behavior of 
the symmetry energy under controlled conditions. To the best of our knowledge, the con¬ 
straints on the snpra-satnration density behavior of the symmetry energ y ob tained so far are 
all from the transport model analyses on the data of 7r“/7r+ ratio |l9l425| and n/p elliptic 
flows 26|, 271 in heavy-ion collisions but unfortunately they are contradictive with each other, 
leaving a confusing situation for the community. 


Conventionally, the nuclear matter EOS is defined as the binding energy per nucleon as 
a function of the density and a number of bulk characteristic parameters defined at satura¬ 
tion density po are usually introduced to quantitatively characterize the energy of symmetric 
nuclear matter and the symmetry energy j28|, [29[. For example, the energy E^o(po) and in¬ 
compressibility Kq of symmetric nuclear matter are the two lowest order bulk parameters 
for the EOS of symmetric nuclear matter while the symmetry energy magnitude ^sym(po) 
and its slope parameter L are the two lowest order bulk parameters of the symmetry en¬ 
ergy. While several lower order bulk characteristic parameters of asymmetric nuclear matter, 
such as Eoipo), ^sym(po) and L have been relatively well constrained or in significant 
progress yet the higher order bulk characteristic parameters are still poorly known. 

Actually, there has little experimental information on the third-order derivative parameter 
Jo of symmetric nuclear matter at po 0 and the symmetry energy curvature parameter 
^sym [29|. However, the higher order bulk characteristic parameters have been shown to 
be closely related to some important issues in nuclear physics and astrophysics such as the 
determination of the isobaric incompressibility of asymmetric nuclear matter [281 . [32| and the 
core-crust transition density and pressure in neutron sta rs [33l 435|. In particular, within the 
Skyrme energy density functional, it has been proposed [29| that the higher-order curvature 
parameter i^sym may play an important role in the determination of the snpra-satnration 
density behaviors of the symmetry energy. 

So far (very likely also in future), essentially all the obtained constraints on E^sym(p) are 
based on some energy density functionals or phenomenological parameterizations of J^sym(p)- 
Therefore, it would be very interesting to see whether there exist some universal laws (system- 
atics) for the density dependence of the symmetry energy within these functionals or param¬ 
eterizations and whether one can get some useful information on the high density symmetry 
energy from the relatively well-known knowledge of the symmetry energy around saturation 
density. For all the energy density functionals or phenomenological parameterizations, the 
^sym(p) increases from p = 0 up to a certain density around po and then either continu¬ 
ously increases or decreases depending on the parameters of the energy density functionals 
or phenomenological parameterizations. While the parameters ^sym(po), L and i^sym accu¬ 
rately characterize the symmetry energy density behaviors around po, their relation to the 
density behaviors at sub- and snpra-satnration densities in various energy density functionals 
or phenomenological parameterizations of E^sym(p) are still unclear. 

In the present talk, we report the preliminary results of the study on the systematics 
of the density dependence of nuclear matter symmetry energy in the ambit of microscopic 
calculations with various energy density functionals. We systematically analyze the relation 
between the parameters ^sym(po )5 L and i^sym defined at saturation density po and the 
symmetry energy density behaviors at sub- and snpra-satnration densities in various energy 
density functionals. In addition, since the dense matter at high baryon densities could be 
quark matter, we also discuss briefly the quark matter symmetry energy. 




International Workshop on Multi facets of Eos and Clustering (IWM-EC 2014) 


2 Symmetry energy systematics 

The EOS of isospin asymmetric nuclear matter, given by its binding energy per nucleon, can 
be expanded to 2nd-order in isospin asymmetry S as 

E{p, S) = Eo{p) + E,y^{p)5‘^ + (1) 


where p = + Pp is the baryon density with pn and pp denoting the neutron and proton 

densities, respectively; 6 = {pn — pp)/p is the isospin asymmetry; Eq{p) = E{p,6 = 0) is the 
binding energy per nucleon in symmetric nuclear matter, and the nuclear symmetry energy 
is expressed as 


-^sym (p) 


1 d^E{p,5) 
2 ! dS^ 


( 2 ) 


Around a reference density pr, the E^sym(p) can be expanded in Xr = {p ~ Pr)/Pr as 


Esym{p) — Esym{Pr) H-H-sym(Pr)^2 O(x^), (3) 

where L{pr) = 3pr^^^§f^\p=p, and K^ym{pr) = ^pl^^f^^\p=Pr are, respectively, the 
slope and curvature parameters of the symmetry energy at pr, and they are the lowest- 
order two bulk parameters characterizing the density behaviors of the symmetry energy 
around pr. In particular, when the reference density pr is taken as the saturation density po, 
the L{pr) and i^sym(Pr) are then reduced to the famous symmetry energy slope parameter 

L = and symmetry energy curvature parameter Xsym = 

respectively. 

To examine the symmetry energy systematics, we select a comprehensive large sample of 60 
well-calibrated interactions in various energy density functionals, namely, 33 Skyrme interac¬ 
tions (v090, MSk7, BSk8, SKP, SKT6, SKX, BSkl7, SGII, SKM*, SLy4, SLy5, MSkA, MSLO, 
SIV, SkSM*, kMP, SKa, Rsigma, Gsigma, SKT4, SV, SkI2, SkI5, BSK18, BSK19, BSK20, 
BSK21, MSLl, SAMi, SV-min, UNEDFO, UNEDFl, TOV-min), 2 Gogny interactions (DIS 
and DIN), 18 nonlinear RMF interactions (FSUGold, PKls24, NL3s25, G2, TMl, NL-SV2, 
NL-SH, NL-RAl, PKl, NL3, NL3*, Gl, NL2, NLl, lU-FSU, BSP, lUFSU*, TMl*), 2 density- 
dependent RMF interactions (DD-MEl and DD-ME2), 3 point-coupling RMF interactions 
(DD-PGl, PC-PKl, PG-Fl), and 2 relativistic HF interactions (PK03 and PKAl). These 
interactions include the 46 interactions used in Ref. jUj (except BCP which is designed for 
density up to only 0.24 fm“^) and other 14 interactions (i.e., BSK18, BSK19, BSK20, BSK21, 
MSLl, SAMi, SV-min, UNEDFO, UNEDFl, TOV-min, lU-FSU, BSP, lUFSU*, TMl*) con¬ 
structed more recently. Shown in Fig. [T]is the symmetry energy as a function of the density 
normalized by the corresponding saturation po with the 60 interactions. It is clearly seen that 
various energy density functionals predict very different density behaviors of the symmetry 
energy, especially at supra-saturation densities. For example, the magnitude of the symmetry 
energy at 2po can be varied from about 15 MeV to 100 MeV, depending on the models and 
interaction parameters. Furthermore, it is seen that some non-relativistic interactions predict 
negative symmetry energy at baryon densities above about 2.5po- 

Around the saturation density po, the magnitude F^sym(p) and the density slope L{p) of 
the symmetry energy can be well approximated, respectively, by 


^sym(p) = ^sym(po) + Ay + 




2 ! 


(4) 
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Figure 1. (Color online) The symmetry energy as a function of the density normalized by the 
corresponding saturation po in various energy density functionals with 60 interactions. See text for 
the details. 


and 

L {p) = Lp/pO + KgyYnXP/ PO: (5) 

with X = {p — pq)/ pQ. Based on a correlation analysis within SHF approach, it has been 
shown [ 2 ^ that ^sym(p) can well describe the magnitude of the symmetry energy up to 2po- 
How well can Eq. Q and Eq. approximate the corresponding values at densities deviated 
from po in various energy density functionals? Shown in Fig.[2]is E^sym(p) vs ^s*ym(p) with the 
60 interactions at p = O.Spo, P = 2^0, P = 2.5po and p = 3po. A very strong linear correlation 

(the Pearson linear correlation coefficient r is larger than 0.98 for all the cases) is observed 

between E^sym(p) and £’sym(p) ^cr the 60 interactions at all the four densities considered here. 
Similarly, Fig. [3] shows the correlation between L{p) and L*{p) with the 60 interactions at 
p = 0.5po, P = 2po, P = 2.5po and p = 3po^ and again a strong linear correlation is observed 
between L{p) and L*{p) (r is larger than 0.93 for the densities considered here). 

The strong linear correlation between E^sym(p) and ^sym(p) well as between L{p) and 
L*{p) for the 60 interactions at different densities shown in Fig. [2] and Fig. [3] leads to the 
following relations 

Esym{p) « A{p) ^ B{p)E*y^{p), (6) 

L{p) « A^{p) + B^{p)L*{p). (7) 

The values of the coefficients A(p), H(p), Al(p) and Hl(p) generally depend the density. In 
particular, one can see from Fig.[2]and Fig.[3]that A{p) (Al(p)) is generally nonzero and B{p) 
(Hl(p)) usually deviates from unit, reflecting the higher-order effects beyond the expansion in 
Eq. Q and Eq. ([5]). We note A{p) ^ 0 (Al(p) ~ 0) and B{p) ^ 1 (Hl(p) ~ 1) for p po as 
expected. Eq. and Eq. © establish the systematics of Esym(p) and L{p) in terms of three 
characteristic parameters E^y^{po)^ L and i^sym- It should be noted that in principe Eq. (O 
can also be deduced from Eq. ^ according to the definition. We would like to point out the 
Eq. dni for the systematics of Esym(p) can be safely applied in the density region from po/5 
to 3po where the Pearson linear correlation coefficient r is always larger than 0.96. Similarly, 
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Figure 2. (Color online) Esym{p) vs E*y^{p) in various energy density functionals with 60 interac¬ 
tions at p = 0.5po, p — 2po, p — 2.5po and p — 3po- 
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Figure 3. (Color online) Same as Fig. but for L{p) vs L*{p). 
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Figure 4. (Color online) Esym(pc) vs E*y^{pc) (a) and L{pc) vs L*{pc) (b) in various energy 
density functionals with the 60 interactions. The constraints Esym(pc) = 26.65±0.20 MeV (AE) and 
L(pc) = 46.0 ± 4.5 MeV (Arnp) obtained in Ref. [sj are also included. 


the Eq. Q for the systematics of L{p) can be safely applied in the density region from po/2 
to 2.5po where the Pearson linear correlation coefficient r is always larger than 0.96. 


3 Supra-saturation density behaviors of the symmetry energy 

The systematics of E^sym(p) and L{p) in Eq. ^ and Eq. Q imply that the three characteristic 
parameters Esym{po)^ L and i^sym (and thus E^sym(p) and L{p)) can be determined once three 
values of either E^sym(p) or L{p) are known. This means that one can extract information 
on the high density behaviors of the symmetry energy from the relatively well constrained 
(sub-)saturation density behaviors of the symmetry energy. 

In recent years, several accurate constraints on the symmetry energy at subsaturation 
density have been obtained through analyzing nuclear structure properties of heavy nuclei. 
Indeed, a quite accurate constraint on the symmetry energy at the subsaturation cross density 
pc = 0.11 fm“^, i.e., Esym{pc) = 26.65 ±0.20 MeV, has been recently obtained from analyzing 
the binding energy difference of heavy isotope pairs js^l- At the same time, an accurate 
constraint on the density slope at pc, ke, L{f^ = 46.0 ± 4.5 MeV has been obtained from 
analyzing the neutron skin data of Sn isotopes [STj. At density pc = 0.11 fm“^, the systematics 
of E^sym(p) and L{p) in Eq. (|6]) and Eq. ([7j) lead to the following expressions 

-^sym(Pc) ~ Oj{Pc) ± ^(Pc)-£^sym(Pc): (8) 

L'^ipc) ~ <^l(Pc) ± ^l(Pc)A(pc). (9) 

The values of coefficients a(pc) and b{pc) (nL(Pc) and 6 l(Pc)) can be obtained from linear 
fitting to the correlation between E^y^{pc) and E^sym(Pc) {L*{pc) and L{pc)). Shown in 
Fig. [4] are E^sym(Pc) vs E^y^{pc) and L{pc) vs L*{pc) with the 60 interactions, and one can 
observe a very strong linear correlation (the Pearson linear correlation coefficient r is about 
0.999 for both cases) between E*y^{pc) and ±^sym(Pc) as well as between L*{pc) and L{pc)^ 
and these linear correlations lead to a(pc) = —0.111 ±0.111 MeV, b{pc) = 1.007 ± 0.004, 
aL(pJ = -4.148 ± 0.358 MeV and 6 l(Pc) = 1-048 ± 0.006. 
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Figure 5. (Color online) Density dependence of Esym{p) and L{p) according to the systematics in 
Eq. dSl) and Eq. 0 with Esym(po) = 32.5 ± 0.5 MeV, L = 46.7 ± 12.8 MeV and Ksym = -166.9 ± 
168.3 MeV. The results from the MDI interaction and the phenomenological parameterizations of 
Esym,pot(p) ~ (p/po)^ are also include for comparison. 


Besides E^sym(Pc) and I/(pc), one needs another constraint condition to determine 
^sym(po), L and i^sym and thus the supra-saturation density behaviors of the symmetry 
energy. In the present work, we further use the constraint of E^sym(po) = 32.5 ± 0.5 MeV 
obtained recently from a new and more accurate finite-range droplet model analysis of the nu¬ 
clear mass of the 2003 Atomic Mass Evaluation j^. Therefore, from E^sym(Pc) = 26.65 ±0.20 
MeV, L{pc) = 46.0 ±4.5 MeV and T^sym(po) = 32.5 ±0.5 MeV, one can obtain L = 46.7 ±12.8 
MeV and iCgym = —166.9±168.3 MeV. It is interesting to see that the obtained L = 46.7±12.8 
MeV is in very good agreement with other constraints extracted from terrestrial experiments, 
astrophysical observations, and theoretical calculations with controlled uncertainties (iMi. 
The obtained i^sym = —166.9± 168.3 MeV also agrees well with the result i^sym = —100±165 
MeV [i^ obtained from a correlation analysis within SHF approach. 

Based on E^ymipo) = 32.5±0.5 MeV, L = 46.7±12.8 MeV and K^ym = -166.9±168.3 MeV, 
one then can obtain E^sym(p) and L{p) according to the systematics in Eq. 0 and Eq. 0, 
and the results are shown in Fig. [5l For comparison, we also include the results from the 
MDI interaction with x = 1, 0, — 1 and —2 and the phenomenological parameterizations 
of ^sym,pot(p) ^ {p/po)^ (jdl. for the potential energy part of the symmetry energy with 
7 = 0.5, 0.75, 1.0, 1.25 and 1.5. We would like to point out that the MDI interaction and 
the parameterizations of ^sym,pot(p) ^ {p/poV have been extensively applied in transport 
model simulations of heavy ion collisions. One can see from Fig. [5] that the present analysis 
based on the symmetry energy systematics with E^sym(po) = 32.5 ± 0.5 MeV, I/ = 46.7±12.8 
MeV and i^sym = —166.9±168.3 MeV favors a softer symmetry energy and suggests that the 
symmetry energy cannot be stiffer than a linear density dependence. 

In particular, at the supra-saturation density of 2pQ^ we find E^sym(2po) = 40.2 ± 12.8 MeV 
and L{2po) = 8.9 ± 108.7 MeV. We note that these values are in nice agreement with the 
variational many-body theory calculation with WFFl interaction (4l| which can give a good 
description on the recent observation of heavy neutron stars with radius of km |42| . 
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4 Quark matter symmetry energy 

At extremely high baryon density, the main degree of freedom could be the deconfined quark 
matter rather than the confined baryon matter, and there the quark matter symmetry energy 
should be involved for the properties of isospin asymmetric quark matter (isospin symmetry 
is still satisfied in quark matter). The isospin asymmetric quark matter could be produced 
in ultra-relativistic heavy ion collisions induced by neutron-rich nuclei and it could also exist 
in compact stars such as neutron stars or quark stars. Although significant progress has 
been made in understanding the density dependence of the nuclear matter symmetry energy, 
there has little information on the density dependence of the quark matter symmetry energy. 
Theoretically, it is difficult to calculate the quark matter symmetry energy since the ab initio 
Lattice QCD simulations does not work at finite baryon density while perturbative QCD only 
works at extremely high baryon density. 

Similarly as in the case of nuclear matter, the EOS of quark matter consisting of u, d, 
and s quarks, defined by its binding energy per baryon number, can be expanded in isospin 
asymmetry 6q as 


where Eo{nB^Tis) = E{nB,Sq = 0, n^) is the binding energy per baryon number in three-flavor 
u-d-s quark matter with equal fraction of u and d quarks; the quark matter symmetry energy 
Esymi'^B^'^s) is expressed as 


A'sym ('^5 5 — 


1 d‘^E{nB,S,ns) 


2 ! 




( 11 ) 


Sa=0 


The isospin asymmetry of quark matter is defined as 


^0 = 3 


Tld Tlu 
nd + Uu 


( 12 ) 


which equals to —nzjnB with the isospin density = Uu — Ud and ub = {riu + ^d)/3 
for two-flavor u-d quark matter. We note that the above definition of 6q for quark matter 
has been extensively used in the literature 43-[4^. and one has 6q = 1 (—1) for quark 
matter converted by pure neutron (proton) matter according to the nucleon constituent quark 
structure, consistent with the conventional definition for nuclear matter, namely, = 

—nsjuB^ In Eq. ([TQ|) . the absence of odd-order terms in Sq is due to the exchange symmetry 
between u and d quarks in quark matter when one neglects the Coulomb interaction among 
quarks. The higher-order coefficients in 5q are shown to be very small in various model 
calculations (t^ . 

It has been demonstrated recently that the isovector properties of quark matter may 
play an important role in understanding the properties of strange quark matter and quark 
stars. If the recently discovered heavy pulsars PSR J1614-2230 Q and PSR J0348+0432 j4§| 
with mass around 2Mq were quark stars, they can put important constraint on the isovec¬ 
tor properties of quark matter, especially the quark matte r sy mmetry energy. Within the 
confined-isospin-density-dependent-mass (CIDDM) model (^, in particular, it has been 
shown that the two-flavor u-d quark matter symmetry energy should be at least about twice 
that of a free quark gas or normal quark matter within conventional NJL model in order to 
describe the PSR J1614-2230 and PSR J0348+0432 as quark stars. 
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5 Summary 

We have explored the systematics of the density dependence of the symmetry energy in the 
ambit of microscopic calculations with various energy density functionals. Our results indicate 
that the symmetry energy magnitude E^sym(p) and its density slope L{p) from subsaturation 
density to supra-saturation density can be essentially determined by three parameters defined 
at saturation density po, be., the magnitude £’sym(po )5 the density slope L and the density cur¬ 
vature ibsym- This finding implies that three values of E^sym(p) or L{p) essentially determine 
Esym{p) and L{p) in large density region. In particular, using E^sym(Pc) = 26.65±0.2 MeV and 
E{pc) = 46.0 ± 4.5 MeV at pc = 0.11 fm“^ extracted from isotope binding energy difference 
and neutron skin of Sn isotopes together with E^sym(po) = 32.5±0.5 MeV obtained from finite- 
range droplet model analysis of nuclear binding energy, we obtain I/ = 46.7±12.8 MeV and 
Ksym = -166.9±168.3MeVas well as E;sym(2po) 40.2 ± 12.8 MeV and L(2po) ^8.9 ±108.7 
MeV. These results favor a soft to roughly linear density dependence of the symmetry energy. 

We have also discussed the quark matter symmetry energy, which has been shown to play 
an important role in understanding the properties of strange quark matter and quark stars. In 
particular, it has been suggested that the two-flavor u-d quark matter symmetry energy should 
be at least about twice that of a free quark gas or normal quark matter within conventional 
NJL model in order to describe the recently discovered heavy pulsars PSR J1614-2230 and 
PSR J0348±0432 with mass around 2Mq as quark stars. 
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hai Municipal Education Gommission and Shanghai Education Development Foundation, the 
Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions 
of Higher Learning, and the Science and Technology Commission of Shanghai Municipality 
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